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We analyze superconducting instabilities in 3D and 2D extended Hubbard model with Coulomb repulsion 
between electrons on neighboring sites in the limit of low electron density (n e ; — > 0) on simple cubic (square) 
lattice. We show that in a realistic strong-coupling case U S> V S> W (U and V are the onsite and the 
intersite Coulomb repulsions, W the bandwidth) the main SC instability corresponds to the p-wave pairing 
and in the leading order is correctly described by the equations obtained earlier in the absence of the intersite 
Coulomb interaction V = 0. 



1. INTRODUCTION 



One of the main challenges of the modern condensed 
matter physics is to identify the origin of superconduc- 
tivity in supcrfiuid 3 He, heavy fermion compounds and 
Si'2Ru04, semimetals and supcrlatices. A lot of the ex- 

Spcrimental data as well as theoretical calculations sug- 
gest that the pairing results from the electron-electron 
interaction. In this scenario, a Coulomb repulsion is in- 
verted into attraction due to the fermion background 
.and retardation effects. This was first suggested by 
Kohn and Luttinger pQ for a 3D system with point- 
like repulsion. The authors of Refs. El H] extended 
^ jthe analysis to 2D systems and took into account the 




o 



effects of long range Coulomb interaction in dense elec- 
tron plasma. Recently, the question about the role of 
■full Coulomb interaction for non-phonon mechanisms 
of superconductivity was raised in connection with the 
HTSC physics by Alcxandrov and Kabanov [5], and it 
still demands very thorough investigations both in the 
jellium and lattice models. 

In the present paper we consider the simplest and 
the most repulsive (thus the most unfavorable for effee- 
■tive attraction and SC) lattice model with the strong 
.on-site Hubbard repulsion U and the relatively strong 
additional Coulomb repulsion V on the neighboring sites 
(Fig. [1]) . We show that in this model the p-wave super- 
conductivity exists in both the 3D and 2D case [2J 02 2] . 
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We assume the following estimates: U 



'/ecis for 



Hubbard U and V ~ e 2 /ed for Coulomb V. Here 
as ~ e/me 2 is the Bohr radius, e the effective dielectric 
permittivity, d the intersite distance. We assume that 
for e ~ 1: a B ~ 0.5 A, and d ~ 3 - 4 A. 

In the simple 3D cubic lattice the bandwidth is 



~d/2 d r 

Fig. 1. Effective interaction in the extended Hubbard 
model with Coulomb interaction on neighboring sites. 

W = 12t where t is the hopping integral, and the 
electron mass at low density (practically, empty lat- 
tice) m = l/2td 2 . The uncorrclatcd electron spectrum 
e(p) = — 2t (cos p x d + cos p y d + cos p z d) approximately 
has quadratic form e(p) = —W/2 + p 2 /2m. Similarly, 
the chemical potential measured from the bottom of the 
band reads (J, = —W/2 + sp, where Sf = p%/2m is the 
Fermi energy, pp the Fermi momentum. If, as usual, 
we assume as <C d (which, rigorously speaking, is valid 
at moderate values of e > 1), then comparing the esti- 
mates for U ~ e 2 /£as, V ~ e 2 /ed, and W ~ 1/md 2 in 
the limit as/d <S^. 1 we come to the following hierarchy 
of parameters: 

U > V > W. (1) 

Note that some important SC systems possibly includ- 
ing HTSC could have large values of e and thus be in 
difficult intermediate regime. 

In this paper we construct the theory for the SC in- 
stability in the parameter range ((T|) and at low electron 
density n e i — !> (or subsequent gas parameter ppd <C 
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1), neglecting an important question of the microscopic 
phase separation of Mott-Hubbard type [6] [7] with FM 
polarons inside the AFM-matrix, and that of Verwey 
type [SI |H] with metallic polarons inside the charge- 
ordered matrix. These instabilities towards nanoscale 
phase separation arise in the model under condition ([1]) 
close to n e i — > 1 for the Mott-Hubbard and n e i — > 1/2 
for the Verwey type of phase separation. 

In the following we show that the leading SC insta- 
bility at n e i —¥ corresponds to the triplet p-wave pair- 
ing and in the leading order of the gas parameter |10j 
p F d is described by the expressions obtained in Rcfs. 
[3] 131 2] for the low density Hubbard model in the ab- 
sence of Coulomb interaction (at 1^ = 0). We review 
the 2D case and present analogous results for the p-wave 
pairing in the strong coupling case [TTJ [TJ] which is also 
in accordance with the low-density Hubbard model in 
the absence of V (at V = 0). 



2. THE MODEL 

We consider the Hamiltonian 
H'=H-fiN = -t c L c j<r + U J2 n * ni ± 
V 



<IJ><7 



(2) 



<IJ> io 

where n.; CT = cj^c^ is the electron density on site i with 
spin projection a. After Fourier transformation, the 
Hamiltonian reads: 

pa pp'q 

+ y (P'P') C P<x C p<+ qCT 'Cp +q a'C p v, ( 3 ) 

pp'qcr a' 

where 

^(P,P') = V[cos(p x -p' x )d+cos(p y -p' y )d+cos(p z ~p' z )d]. 

(4) 

In analogy with Ref. [Tj3] it is useful to expand the ef- 
fective interaction U e ff = U + 2V(p, p') into the sum 
of the s-wave and p-wave partial harmonics. 

At the low density pd <C 1 the expansion up to 
quadratic terms gives effective interactions for s-wave 
and p-wave harmonics correspondingly: 



U, 



eff = U + 6V + o(p 2 d 2 ), and U p eff = 2Vpp'd 2 . (5) 



In the strong-coupling case U 3> V 3> W it is conve- 
nient to renormalize U*fj and in terms of vacuum 
Kanamori T-matrices T s and T p [2]. To do that we 



solve the Bethe-Salpeter equation in vacuum [TS]. This 
yields [4] in the low-energy sector: 

(U + 6V)d 3 (U + 6V)d 3 



i + (u + 6V)d 3 J^^ + 



(6) 



where j3 s ~ ^ U g^^ > is the Born parameter for the s- 
wave channel, and we neglect the antibound state which 
corresponds to the pole of the T-matrix at high energies 

E~ U [SUB]. 

We can introduce the s-wave scattering length 

mT s T s j3 s d 



a s 



4tt 



8ntd 2 (1 + /3 S ) ' 



(7) 



and in the strong-coupling limit j3 s 3> 1, evidently, 
a s ~ d (see Ref. g]). 

Correspondingly, the 3D gas parameter of Galitskii 
[Uj: 

^ _ 2a s p F ^ 2p F d ^ 

7T 7T 

Note that the same result for the s-wave scattering 
length is valid in the strong-coupling low-density limit 
of the Hubbard model without Coulomb interaction 

v = o. 

Similarly, for the T-matrix in the p-wave channel 
T p = 2A p pp'd 2 



we get 
A i) 

where 



Vd 3 



Vd 3 



v 



(i + &)' 



11 w >0 



(9) 
(10) 

(11) 



is the dimcnsionlcss Born parameter for the p-wave 
channel and n ~ 1 is a numerical coefficient. 
Introducing the p-wave scattering length 

K Vd 



p 8irtd 2 8nt(l + p p ) 
we obtain in the strong coupling case (3 p ^ 1 
a p ~ d. 

Thus 



mT, 



4tt 



- = 2a p pp'd 2 ~ d 3 pp' ~ d 3 pp' cos 0, 



(12) 

(13) 
(14) 



where 9 = pp', and thus the dimensionless p-wave gas 
parameter reads: 

p F d 3 p F ~ {p F df. 



(15) 

Note that the estimate (fT5"j) is natural for the p-wave 
harmonics of the scattering amplitude for slow (p F d < 
1) particles in vacuum |15j . 
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Fig. 2. Irreducible diagrams in the second order of the 
s-wave gas parameter \ s which are nonzero only in the 
presence of fermion background (at ef 0). 

5. BETHE-SALPETER INTEGRAL EQUATION 
FOR T c 

According to Landau- Thouless criterion for SC [17] . 



1 + r ; In 



(16) 



where C ps 0.58 is the Euler constant, T the total vertex 
for the Cooper channel, V the irreducible bare vertex, 
and I the orbital moment of the Cooper pair. 

The critical temperature T c is given by the pole of 
(fT!)]) . If f < for several values of I, then the actual 
symmetry of the superconducting state corresponds to 
the highest T c . According to Kohn and Luttinger pQ, in 
the absence of the Coulomb interaction (i.e. of X p ([T5|) ) 
T/^o is given by the sum of four irreducible diagrams 
(see Fig. [2J which are of the second order of the s-wave 
gas parameter A s . 

At the same time, for I = 1 due to the presence of 
Coulomb repulsion V: 



r*=i = A;?nf =1 + x p 



(17) 



where n d = LT(p + k) is the exchange diagram (see Fig. 
(2]d.), X p is the bare vertex due to the p-wave vacuum 
contribution of the intersite Coulomb interaction V. 

As shown in Ref. [T] for contact interaction A s the 
first three diagrams in the Fig. [5] exactly cancel each 
other, and the resulting Ti = i is given by the fourth, ex- 
change diagram (see ([17])). 

An exact evaluation of simple integrals shows 013111] 

A 2 

that for the exchange diagram A^LT^ = — < 
which corresponds to the attraction and cannot be over- 
compensated by the repulsive bare vertex contribution 
Xp r~j X 3 . This contribution only changes the next term 
in the expansion of T; = i in terms of gas parameter and, 






Fig. 3. Irreducible diagrams in the third order of gas 
parameter for the Cooper channel. 

in fact, is the corrections to main exponent. To be spe- 
cific (see Ref. PHE] and Fig. 



Ti= 



As 


(1- 


1_ 13 





X p +o(Xt). 



(18) 



Let us repeat again that the result (jTHJ) with the 
bare vertex A p ~ X 3 ~ (ppd) 3 is to some extent evi- 
dent because according to quantum mechanics [TS] for 
slow particles in vacuum the p-wave harmonic of the 
scattering amplitude is of the order (app) 2l+1 ~ (app) 3 
at I = 1 and p ~ pp. Thus the repulsive term in our 
case does not overcompensate the Kohn-Luttinger at- 
tractive contribution which arises only in fermion sub- 
stance (when sp ^ 0) and is proportional to (app) 2 . 
The only peculiarity of the lattice is that a ~ d at the 
large Hubbard U 3> W and at low density of electrons, 
and hence Ag ~ (ppd) 2 <C 1. 

Thus even at the most repulsive (and thus unfavor- 
able for effective attraction and SC) hierarchy of pa- 
rameters U ^> V ^> W the presence of the Coulomb 
repulsion V does not change the main exponent for 
the p-wavc critical temperature which reads T c % ~ 
e^exp (-j?) as in [3 EE]. 

Note that if we change the hierarchy of parameters 
and make Coulomb repulsion weaker W S> U 3> V, 
then in the Born case: 



An 



mVdr 

47T 



3 j3 



-p F d 



mVd 2 . 3 V , 2 
— A^ - — A^ < X 

47T W 



U + 6V 
W 



(19) 

for pp — > and still the overcompensation of the Kohn- 
Luttinger attraction by the bare repulsion due to the 
intersite Coulomb interaction V is impossible. Thus, 
in the principal approximation in the gas parameter we 
restore the results on the possibility of the p-wave su- 
perconductivity obtained earlier in the absence of the 
intersite repulsion [2] [3] 0] . 



4. 2D EXTENDED HUBBARD MODEL 

In the 2D extended Hubbard model with attractive 
interaction (—V < 0) on neighboring sites the vac- 
uum T-matrices for the s-wave and p-wave channels 
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were obtained in the Refs. [T3J E5]. After the sub- 
stitution — V — >• V they yield for the s-wave channel 
U*jf = U + AV in the repulsive case U > V > W: 



An 



l + ([7 + 4F)(^)ln(^) 



(20) 



where 7 ~ 1 is the numerical coefficient. Again, we 
assume that we are in the low energy sector when one 
can neglect the second pole of the T s which corresponds 
to the antibound state E ~ U US]- In the Eq. 
(|2"0"j) 14 7 = 8t for the 2D square lattice, and the energy 
E = E + W is measured from the bottom of the band. 
If (U + AV)/W < 1, then 



mT s (E) 



1 



An 



ln (w) 



(21) 



In the Cooper problem \E\ = 2e F and with the loga- 
rithmic accuracy we restore the 2D dimensionless gas 
parameter of Bloom [50] : 



f. 



1 



(22) 



where n = p 2 F /2n is the electron density in 2D. 

Analogously, in the p-wave channel U%ff 
2Vpp'd 2 , and the p-wave T- matrix reads: 



mT„ 2mA. 



An 



An 



pp'd 2 



where 



pp', and 

mA p 



2mA f 
An 



iVd 2 



pp'd 2 cos< 



Air (1 + V/V cp )8n 
Correspondingly [TO] . 



V„ 



11.2* « 1.4W. 



(23) 



(24) 



(25) 



At F ^> V^p, the dimensionless p-wave scattering length 
in 2D reads 

mAr, 



p 



■ Cp 



An 



8 n 



(26) 



and, accordingly, the dimensionless p-wave gas parame- 
ter is 

U-^^-pU'-pU 2 - (27) 

Thus f p ^ p 2 pd 2 again in agreement with general 
quantum-mechanical results [15] for slow {p F d < 1) par- 
ticles in vacuum in the 2D case. 



5. THE COOPER PROBLEM IN 2D AT LOW 
ELECTRON DENSITY AND IN THE 
PRESENCE OF INTERSITE COULOMB 
REPULSION 

If we restrict ourselves to a very low electron den- 
sity n e \d 2 <C 1 and quadratic spectrum e(p) — p = 
(p 2 — p 2 F )/2m, then in the second order of the s-wave gas 
parameter the irreducible vertex for the Cooper channel 
reads: 

f = / s 2 n( P + k) (28) 

However, the specific form of the polarization operator 
on quadratic spectrum in 2D [21] for q = p + k 



n(g) = 1-ReWl 



4P 2 , 



(29) 



Px+Py) d 2 /2Am from the quadratic 
+ p 2 /m, sec Rcf. [22]). At very low 



makes the large Kohn's anomaly ineffective for the SC 
problem [1J[TTJ. Indeed, in the SC problem q < 2pp, 

Re^/ 1 — = 0, and thus 11(g) = 1. Hence, the po- 
larization operator does not depend on q, and corre- 
spondingly it does not contain harmonics with / ^ 
(or more precisely, with the magnetic quantum num- 
ber m 7^ 0). Thus n m= x = 0, and SC arises only 
in the third order of f s for quadratic spectrum (or in 
the second order of f s if we take into account correc- 
tions (px+Py)d 2 /m which differ the exact spectrum 
on the square lattice e(p) = — 2t(cosp x d + cosp y d) 
-W/2 +p 2 /m - 1 
one e(p) = — W/2 
density n e i — > the third order terms in the quadratic 
spectrum from three irreducible diagrams in the Fig. [3] 
dominate over the quartic corrections to the spectrum. 

Chubukov [11] found the leading contribution to 
r m= i from the first skeleton diagram in which the 
Cooper loop is inserted into the polarization loop (it 
is important that this diagram is still irreducible with 
respect to Cooper channel). Moreover, the character of 
the large 2D Kohn's anomaly in this diagram changes 
and it becomes Key / 2p F — q. Thus, the Kohn's anomaly 
becomes effective for SC in the third order. As a result 
he has obtained r m =i = —4.1/^ in the Ref. [TTj. In 
the Ref. [12], all three irreducible skeleton diagrams on 
Fig. [3] were calculated numerically on equal ground. As 
a result, the exact vertex 



r m=1 = -6.i/ s 3 



(30) 



is even a little bit more attractive. The details of this 
calculation will be published in a separate article. 
Thus, the total T m= i at n e i — > reads: 



-6.1ft + ap F d 2 + o(ff) 



(31) 
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where a ~ 1 is a numerical coefficient. 

Of course, keeping in mind that f s ~ 1/ In (l/nci 2 ) ~ 
f/ln (l/p^d 2 ), we see that / s 3 > p|d 2 at p F d < f. 
Thus, r m= i « —6.1/1 just like in the case V = 0. 

We can see again that in the strong-coupling limit 
U V ^s> W of the extended Hubbard model on the 
square lattice and at low electron density ppd < 1 an in- 
clusion of Coulomb repulsion does not change the main 
exponent for the p-wave critical temperature 

T cl ~ £F exp(-^). (32) 

Thus in the principal order in the gas parameter we 
again restore the results on the p-wave superconductiv- 
ity obtained earlier [TTJ[T2] in the absence of the intersite 
Coulomb repulsion. 



6. DISCUSSIONS: THE CASE OF LARGER 
DENSITIES 

If we increase the density in the 2D case, we should 
remember that at U V ^> W the homogeneous 
metallic state stretches only up to the density n e i = 
\ - S c , where in 2D S c ~ (W/V) 1/2 (see Ref. [9]). At 
n e i > h — S c the system undergoes a phase transition 
into phase-separated state with metallic clusters inside 
charge-ordered checkerboard matrix (see Fig. 2]). 

Note that at n e i = 1/2 (quarter- filled band) we have 
Verwey localization (charge ordering) due to the con- 
dition V W. Thus, we cannot extend our calcu- 
lations for T c in homogeneous case to densities larger 
than n e i = 1/2. However, it is interesting to construct 
the SC phase diagram of the extended Hubbard model 
with the Coulomb repulsion on neighboring sites at the 
intermediate density n e i < 1/2, and to find the regions 
that correspond to the p-wave, d xy , and d x i_ y 2 pairings 

[31221 123 • 

Another interesting question would be to add to 
the model an infinite set of Coulomb repulsion terms 
with the amplitude decreasing with the distance be- 
tween the sites: V2Tiini+2 on next-to-nearest sites with 
V% < V , VsniHi^ on next-to-next-to-nearest sites with 
V3 < V~2 < V etc. and to build a bridge between 
the extended Hubbard model and the jellium model for 
screened Coulomb interaction considered in Ref. [5J. 

We think, however, that at least at very low elec- 
tron density n e i — > our results on the p-wave critical 
temperature will be stable in the main order of the gas 
parameter ppd -C 1 in 3D and 1 / \n(\ / p F d 2 ) in 2D. 



(pOxOQ)xOxOxO(pOxOd) 

Fig. 4. Phase separation at the density n e ; < 1/2 into 
metallic droplets in charge-ordered matrix. 

7. CONCLUSION 

We considered the extended Hubbard model with 
Coulomb repulsion on the neighboring sites in the most 
repulsive (and thus the most unfavorable for effective 
attraction and SC) strong-coupling case U 3> V ^> W. 
In the limit of small electron density ppd -C 1 we found 
that the contribution from the intersite Coulomb repul- 
sion V to the irreducible bare vertex T/ = i in the p-wave 
channel is proportional to (ppd) 3 in 3D and to (ppd) 2 
in 2D in agreement with general quantum-mechanical 
results for slow particles in vacuum. 

Thus both in 3D and 2D these repulsive terms can- 
not ovcrcompcnsatc attractive contributions which are 
proportional to (ppd) 2 in 3D and to 1/ ln 3 (l/p|nd 2 ) in 
2D. Note that the attractive contributions appear only 
in the presence of fcrmion background (ep ^ 0). Thus 
the results of Refs. [1 QU Q2 [T£] on the p-wave SC of 
Kohn-Luttinger type [S] both in 3D and 2D repulsive- U 
Hubbard model at low electron density and strong cou- 
pling U ^> W are robust against the addition of even 
strong Coulomb repulsion on neighboring sites V ^> W 
in the extended lattice models. Hence we can see that 
the p-wave superconductivity exists in purely repulsive 
models without electron-phonon interaction. 

Note that we can strongly increase the p-wave crit- 
ical temperature already at low density in a spin- 
polarized case [21] or in the two-band situation [35] 
and thus reach the realistic values of T c (of the or- 
der of 1 — 5 K especially in the 2D or in layered sys- 
tems [26j). The p-wave pairing is realized or can be 
expected in superfluid 3 He and ultracold Fermi-gasses, 
heavy fermion compounds and Sr2Ru04, semimetals 
and supcrlatices, layered dichalcogcnides and organic 
superconductors [271 [25] [2^] |5D] . 
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